arXiv:1506.03982vl [math.AP] 12Jun2015 


ON SOME NONLINEAR FRACTIONAL EQUATIONS INVOLVING THE 

BESSEL OPERATOR 

SIMONE SECCHI 
Dedicated to Francesca 


Abstract. Under different assumptions on the potential functions h and c, we study the fractional 
equation {I — A)°‘u = \h{x)\uY~^u + c{x)\u\‘^~^u in R^. Our existence results are based on 
compact embedding properties for weighted spaces. 


1. Introduction 

In this paper we provide some existence results for a class of nonlinear fractional equations of 
the form 

(1.1) (/— A)“ tt = A6(x)|u|^“^u + c(x)|u|'^“^u in 

where 0 < a < 1, p and q belong to the interval (0, 2*) with 2* = 2N/{N — 2a). We will assume 
throughout that 

(H): The potential functions b and c are continuous and bounded. 

Equations involving the nonlocal operator (/—A)“ arise in the study of standing waves ijt = x) 
for Schrodinger-Klein-Gordon equations of the form 

= (^ - (Lx)eMxM^ 

which describe the behavior of bosons. We refer to [22, 23] for a physical introduction to these 
fractional equations. 

The case a = 1 corresponds to the classical Schrodinger equation —An + n = f{x,u), and we 
cannot review the huge literature here. In the last years the so-called fractional Laplacian (—A)" 
has become popular in the community of Nonlinear Analysts: we refer the interested reader to 
the recent survey [6] and to the references therein. 

The most physically relevant fractional case is a = 1/2, and the corresponding Bessel operator, 
or more precisely the operator \/—A + — m, goes under the name of relativistic Schrodinger 

operator. 

It is interesting to quote a sentence from [7, page 119]: 

One of the fundamental mathematical problems in proving the stability (or in¬ 
stability) of matter is to estimate the infimum of the spectrum of the operator 
H = Hq -|- V when the numbers N and/or M become large. In this asymptotic 
regime, the free Hamiltonian Hq = F{p) = V—A -|- — m can be abandoned 
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and replaced by Hq = |p|. Indeed, the difference between these two operators re¬ 
mains bounded and the asymptotic result needed to prove the stability of matter 
can be proved using either one of these free Hamiltonians. Obviously, the scaling 
properties of the function F{p) = |p| attracted early investigations of the corre¬ 
sponding pseudo-differential operator Hq especially because its scaling is related 
to the scaling of the Coulomb potential. Also, very fine estimates on its Green’s 
function are available. These are the technical reasons why F[p) = \p\ is preferred 
to \/p^ + w? — m. 

We will come back to the issue of scaling properties in the last Section. R. Frank et al. in 
[18, page 5] promise to extend part of their spectral theory for the fractional Laplacian to the 
pseudodifferential operator (—A -|- m^)". 

The equation y/I — /S. u = f{u) was studied in [37] by means of a Dirichlet-to-Neumann local 
realization that was extended to any a S (0,1) by Fall et al. in [15]. The case a = 1/2 was also 
studied in [11-14] with a non-local convolutions term on the right-hand side. 

On the contrary, very few papers deal with equation (1.1) for arbitrary a. € (0,1). The very 
recent paper [17] deals with the case in which the nonlinearity / = f(x,u) is essentially of the 
form f{u) + a(a:)(]M] -|- |u|^) with lim| 2 ;|_>,_|_ooa(ic) = 0, in the spirit of [30]. 

In this paper we limit ourselves to a somehow particular class of nonlinearities as in (1.1), and 
we provide a few existence results under different assumptions on the two potential functions b 
and c. 

Section 2 contains some preliminaries on the functional setting that we use to solve (1.1). In 
Section 3 we begin with the case p>2, q>2,b = 0 identically and c positive and “vanishing at 
infinity” in an appropriate sense. In Section 4 we allow b to tend to a positive constant at infinity, 
and in Section 5 we consider the concave-convex case 1 < p < 2 < q with sign-changing potentials 
b and c. 


Notation. 

(1) The letter C will stand for a generic positive constant that may vary from line to line. 

(2) The symbol || • ||p will be reserved for the norm in LP(M.^). 

(3) The operator D will be reserved for the (Frechet) derivative, also for functions of a single 
real variable. 

(4) The symbol will be reserved for the Lebesgue A^-dimensional measure. 

(5) The symbol CA will denote the complement of the subset A (usually in R-^). 

(6) = {(x,y) SR^ X [0,+oo)} 

(7) For a real-valued function /, we set = max{/, 0}, the positive part of /. The negative 
part of / is dehned similarly. 

(8) The Fourier transform of a function / will be denoted by Fu. 


2. Preliminaries and functional setting 
For a > 0 we introduce the Bessel function space 

^ 0,2 (rAT) = |/: f = Ga * g ioT some g E L^(R^)| 
where the Bessel convolution kernel is defined by 


( 2 . 1 ) 


Ga{x) = 


1 


(47r)“/2r(a/2) Jo 


TT 

exp [ 


exp(--)t 2 


^ dt 
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The norm of this Bessel space is ||/|| = ||( 7||2 if / = Ga * g- The operator (/ — A)~°‘u = G2a *u is 
usually called Bessel operator of order a. 

In Fourier variables the same operator reads 


( 2 . 2 ) 

so that 


= 


J-lo(^(l + |^|2) 


-a/2 






For more detailed information, see [1,33] and the references therein. 


Remark 2.1. In the paper [15] the pointwise formula 

(2.3) {I - A)°‘u{x) = CN,aPA. [ — ^^^K n+2c (\x -y\)dy + u(x) 

was derived for functions u S C'^(M^). Here C 7 v,a is a positive constant depending only on N and 
a, P.V. denotes the principal value of the singular integral, and is the modified Bessel function 
of the second kind with order u (see [15, Remark 7.3] for more details). Since a closed formula for 
Ki, is unknown, equation (2.3) is not particularly useful for our purposes. 


We recall the embedding properties of Bessel spaces (see [16,33,35]). 

Theorem 2.2. (1) = iy“’2(M^) = H“(M^). 

(2) If a > 0 and 2 < (/ < 2* = 2N/{N — 2a), then L"’^(R'^) is continuously embedded into 

if2<q<2*^ then the embedding is locally compact. 

(3) Assume that 0 < a < 2 and a > N/2. If a — N/2 > 1 and 0 < g < a — N/2 — 1, then 
L“’2(E'^) is continuously embedded into C^’^(M'^). If a — N/2 < 1 and 0 < g < a — N/2, 
then L“’^(M'^) is continuously embedded into 

Remark 2.3. Although the Bessel space L“’^(M^) is topologically undistinguishable from the 
Sobolev fractional space we will not confuse them, since our equation involves the 

Bessel norm. 


We collect here a couple of technical lemmas taken from the remarkable paper [29]. 

Lemma 2.4. Let 0 < a < N/2 and let u G L“’^(M^). Let ip G and for each R > 0 

let p>r{x) = (p{R~^x). Then Huir^q (pRU = 0 in L“’^(M'^). If, in addition, ip equals one in a 
neighborhood of the origin, then lim/j_>.+oo Pru = 0 in L“’^(R'^). 

Proof. Since L"’^(R'^) is equivalent to H®(R'^), the proof of [29, Lemma 5] carries over with only 
minor modifications. □ 

Lemma 2.5. Let 0 < a < N/2, let LI be a bounded open subset of and let p G C'“(R'^). 
Then the commutator [p, [I — A)“/^]: L"’^(R'^) L^(R'^) is a compact operator. 

Proof. It suffices to remark that the proof of [29, Lemma 6] actually contains a proof of our 
statement. □ 


Definition 2.6. We say that u G L“’^(R'^) is a weak solution to (1.1) if 

/ (/— A)“/^u (I — A)"/^udx = / b{x)\u\^~‘^uv dx + [ c{x)\u\'^~‘^uv dx 

JRN Jrjv Jrjv 

for all V G L"’^(R'^), or, equivalently, 

/ {1 + \^\‘^)°‘iFu{^)iFv{^) d^ = [ b{x)\u\^~'^uv dx + [ c{x)\u\'^~'^uv dx 

JRN Jrjv J^n 
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It is readily seen that weak solutions to (1.1) correspond to critical points of the differentiable 
functional J: L"’^(]R'^) —)■ R defined by 

J(ri) =11(7 — A)“/^n ||2 — — / b{x)\u\^ dx — ^ [ c{x)\u\‘^ dx 

2 p Jrn q Jrjv 

3. Potentials vanishing at infinity 

In this section with treat the particular case of equation (1.1) with 6 = 0, or A = 0. We can 
actually deal with a more general model, 

(3.1) (7 — A)"tt = c(x)/(tt) in R^ 

Definition 3.1. We say that a function c S L°°(R'^) satisfies the compactness condition (K) if 

lim / c{x) dx = 0 

7A„nCB{0,r) 

uniformly with respect to n whenever {An}n is a sequence of Borel sets such that 

sup£'^(74„) < oo. 

n 

Following [2], we assume 

(cl): c > 0 everywhere. 

(c2): c satisfies the condition (K). 

Our assumptions on the nonlinearity / read as follows: / is continuous and 
(fl): limsup, .n/(s)/s = 0. 

(f2): limsup^^+oc/(s)/s2--i = 0. 

(f3): The map s i-A s~^f{s) is increasing and the primitive F of f satisfies lims_>.+cxD T('S)/s^ = 
+ 00 . 

Remark 3.2. Since Ga > 0, it follows easily that (7—A)“ satisfies a maximum principle. Therefore, 
we will look for positive solutions and assume that f{s) = 0 for every s < 0. 

The Euler functional associated to (3.1) is J: L“’^(R'^) —>■ R, defined by 

= -I c{x)F{u) dx. 

2 Jrjv 

The following lemma can be proved by standard techniques. 

Lemma 3.3. The functional J is of class has the mountain-pass geometry. 

As a consequence, there exists (see [8]) a Cerami sequence {un}n C L"’^(R^) for J, i.e. a 
sequence such that 

(3.2) lim J{un) = c and lim (1 + \\un\\)DJ{un) = 0, 

n—)-+co n—>-+00 

where 

Cmp = inf sup J{-f{t)) 

7Gr o<t<i 

and 

T = {7 G C([0, 1 ], L ^ V ) I 7(0) = 0, 7(7(1)) < 0} . 

In view of Remark 3.2, we may assume that Un >0 for every n G N. 
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Definition 3.4. For a positive and measurable function w: 


pN 


[0, +oo) we denote 


Proposition 3.5. The space is compactly embedded into L'^{cdC^) for every 2 < q < 

9* 

‘ 

Proof. We fix (7 E (2, 2*) and pick e > 0. There exist sq < si such that 

c(x)|sl'? < eC + |sp“) + Cc(x)x[so,si](l'S|)|'S|^“ 
for every s E M. Integrating, 


/CS( 0 ,r) 


c(x)|tt(x)|'^ dx < eCQ{u) + C 


AnCs(o,r) 


c{x) dx 


for every u E where Q{u) = \\u \\2 + \\u\\^ 2 * ^ = |x E 

If Un ^ 0 weakly in then 


pN 


So < I'«(x)| < Sij. 


(/— < Ml and ||u„|| 2 ? < Mi 


for some Mi > 0. In particular, the sequence {Q{vn)}n is bounded in M. On the other hand, if 
we set 


then 


An = {a; € I So < |Pn(a:)| < sij , 

So“T^(An) < [ \vn\'^°‘dx < Ml 
J An 


for all n E N. This shows that sup„T^(An) < oo and from our assumption on c there exists a 
number r > 0 such that 

c{x) dx < 


' AnnC-B(0,r) 


We are ready to conclude that 


/ c(x)|un(x)|'^ dx < eCMi + Si“ 

JB(0,r) 

for all n E N. From the compact embedding of into and the boundedness of c, 

we know that 


/AnnCi?(0,r) 


c(x) dx < {CMi + l)e 


Hence Xn 0 strongly in L^{cdC^). The proof is complete 


lim 

n^ + CxD ^ 

N\ 


B{0,r) 


c(x)|xn(x)|'^ dx = 0. 


□ 


Remark 3.6. If assumption (c2) is replaced by lim|a,|_j,+oo c(x) = 0, then the compactness of the 
embedding is classical: see for instance [36, Lemma 3.2]. A different assumption ensuring the 
compactness of the weighted embedding appears in [9]. 

Proposition 3.7. Ifvn^v weakly in then 

(3.3) lim / c{x)F{vn{x)) dx = [ c{x)F{v{x)) dx 

n^+oo J^N J^N 

and 

(3.4) 


lim / c{x) f {Vn{x))Vn{x) dx = / c{x)f{v{x))v{x)dx 

n^+co 7 rjv J’^n 
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Proof. Let us fix g E (2, 2*) and £ > 0. From (fl)~(f3) there exists a constant C > 0 such that 
(3.5) |c(x)F(s)| < eC + c{x)\s\'^ 


for every s E M and every x E 


pN 


lim 

n^+oo 


Proposition 3.5 tells us that 
c{x)\Vn{x)\'^ dx = / c{x)\v{x)\'^ dx, 


so that for some r > 0, 

/ c{x)\vn{x)\'^ dx < £ for all n E N. 

./CB(0,r) 

But {vn}n is bounded in and therefore also in and in We choose 

M 2 > 0 such that /jjat Vnix)"^ dx < M 2 and /j^iv \vn{x)\'^°‘ dx < M 2 for all n E N. Integrating (3.5) 
we get 


/CB(0,r) 


c{x)F{Vn{x)) dx 


< ( 2 CM 2 + 1) £ for all n E N. 


Applying a general convergence result (see for instance [10, Lemma 2.4]) to the sequence , 

together with our assumption (f2) we conclude that 


lim 

n^ + OD 


S(0,r) 


c{x)F{Vn{x)) dx = 


/S(0,r) 


c{x)F{v{x)) dx. 


This shows the validity of (3.3), and a similar argument proves also (3.4). 
Proposition 3.8. The sequence {un}n introduced in (3.2) is bounded in L“’^(M'^). 


□ 


Proof. For every n E N, let E [0,1] be chosen so that J{tnUn) = maxo<t<i J(tUn)- Let us 
prove that the sequence {J{tnUn)'\n is bounded from above in M. The conclusion is trivial if either 
= 0 or = 1. If 0 < < 1, then DJ{tnUn)un = 0. As a consequence, 

‘2J{fn'^ri) — ‘2J(t.}iUY() D J{tYiUjfjUn — / ci^X^Hdx., 

Jmn 

for H{s) = sf{s) — 2F{s). Since Un > 0 and H is non-decreasing, 

‘2J(tn'^n) — I c{x^F[^Un(x^) dx — 2J{Un) DJ(^Un)Un — 2T(u^) T o(l). 

This shows that {J{tnUn)}n is bounded from above. To complete the proof, we argue by contra¬ 
diction. Let us assume that (possibily along a subsequence) lim„_).+cxD II(-^ ~ 2 \)"/^u „||2 = -|-oo. 
We normalize Un by introducing Wn = Un/\\{I — A)“/^nn|| 2 . Without loss of generality, we may 
assume that Wn w weakly in L“’^(M^). Divinding out the relation J{un) c, we find 


0(1) + ^ = / 

2 7® 


c{x)F{Unix)) 

1 ( 7 - A )“/2 


■ dx = 


'«ll2 


f c{x)F{Un{x)) 

IrN |Un(x)|2 


Wn{x)'^ dx. 


By (f3), to each T > 0 we can attach ^ > 0 with the property that |s| > implies F{s) > Ts^. 
Therefore 

1 , r c{x)F{Un{x)) 


o(l) + - > f 

application of 
letting T —)• -|-oo. 


^-Wn{xfdx>T [ c{x)wnixfdx. 

i( 0 )n{|«„|>^} \un{x)\^ 7c«j-i(0)n{|«„|>a 

An application of Fatou’s lemma yields ^ > T /[;^-i(o) c{x)w{x)‘^ dx and consequently tc = 0 by 
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Finally, given T > 0, we have 

J{tnUn) > j ( liTT-= J{TWn) = ^ / c{x)F{TWn{x)) dx. 

\\\{I - A)°‘/^Un\\2 J ^ 

Proposition 3.7 yields lim„_>.+oo /rjv c{x)F{Twn{x)) dx = 0 and thus 

T2 

liminf . 

n^+oo Z 

Letting T —>■ +oo we contradict the boundedness from above of the sequence {J{tnUn)}n- The 
proof is complete. □ 


Let us state and prove the main result of this section. 


Theorem 3.9. Suppose that (fl)-(f3) and (cl)-(c2) hold. Then equation (3.1) possesses at least 
one positive solution. Moreover u G for some qo G [2, +oo) and p, G (0,1). 


Proof. We first prove that our Cerami sequence {un}n is relatively compact in 
aim, by virtue of Proposition 3.8, we may assume that, up to a subsequence, Un - 
Since {un}n is a Cerami sequence, we find 


lim 

n—>-+oo 


- A1“/2 


{I-A) 


= lim 

2 n^+oo ^ 


c{x)f{Unix))Unix) dx. 


By Proposition 3.7, 


To this 
u weakly in 


lim 

n^+oo 


c{x)f{Un{x))Un{x)dx 


c{x) f {u{x))u{x) dx. 


Exploiting again the fact that DJ{un)u ^ 0, we have 

/o 2 

(/ - A)“/2 


U 


c{x) f {u{x))u{x) dx. 


In particular lim„^+oo ||(L - A)"/2^^||2 = \\{I - A)"/2^||2^ OT Un ^ u strongly in But 

J{u) = Cmp and DJ{u) = 0, so that u weakly solves (3.1). The positivity of u follows from the fact 
that Un > 0 and the positivity of the Bessel function Ga see [17, Proposition 3.2]. The regularity 
of u follows with minor changes from the arguments developed in [16]. □ 


Remark 3.10. We have been sketchy about the regularity theory for our solutions, since the Bessel 
operator is precisely the main tool to develop a regularity theory for the fractional Laplacian, 
see [16, Appendix A]. In this sense, the fractional Laplacian is harder to analyze. Moreover, by 
similar arguments as those in [17] it can be shown that our solution decays exponentially fast at 
infinity. This is a common feature for local elliptic partial differential operators, while it is false 
for the fractional Laplacian, see [16]. 


4. Potentials having a finite limit 

The main tool that we used to solve equation (3.1) is the compactness of the embedding 
L"’^(M'^) L^{cdC^) stated in Proposition 3.5. In this section we study a model case in which 
a different approach must be used. We consider 

(4.1) {I — A)°'u = \h{x)\uf’~‘^u + c{x)\u\^~‘^u in M'^, 

where p, q G (2, 2*) and A > 0 is a parameter. We will assume that c > 0 satisfies the condition 
(K), and impose the following ones on the potential 6: 

(bl) b G b > 0 but not identically zero, and lim|a.|_^+oo b{x) = b. 











SIMONE SECCHI 


Weak solutions to equation (4.1) correspond to critical points of the Euler functional J\: 
M defined by 




— — [ b(x)\u(xW dx — - [ c(x)\u(xW dx. 

2 p Jrn ^ q J^N 


Remark 4.1. The functional f c(x)|n(x)|'^ dx is weakly sequentially continuous by the results 
of the previous section. 

Let us introduce the artificial constraint 

Mx = [ue \ {0} I DJx{u)u = O} , 

and standard arguments show that is a natural constraint for Jx. In particular, any solution 
of the minimization problem 

= inf J\{u) 

u^Mx 

is a solution to equation (4.1). 

For A = 0, we introduce the constraint 

Mo = |n s L"’2 (]r^) \ { 0 } I DJo{u)u = o} 

corresponding to the Euler functional 

1 




q JRN 


c(x)|ti(x)|'^ dx. 


By Remark 4.1, the minimization problem 


Jo = inf Jo(ti) 

u£A4.o 


is solved by some function uq S L“’^(R'^) that satisfies {I — A)"uo = c(x)|tto|'^“^tto- 
Lemma 4.2. There results Ix < Iq for all A > 0. 

Proof. Let u G L"’^(R^) be such that 

(/— A)"/^n = [ c{x)\u{x)\^ dx. 

2 JrN 

Pick a G (0,1) such that v = du G M.x- If we differentiate 

h(a) = — (/— A)"/^tt — — / c{x)\u{x)\'^ dx 

2 2 g Jrat 

and remark that Dh{a) > 0 for every a G (0,1), we may conclude that 


h{v) = - (/-A)“/2u 


nP 


a- 


p Jrn 


A6(x)|u(x)|^ dx — 




q Jrn 


c(x)|u(x)|'^ dx 


a 

2 


(/— A)"/^m — — / c{x)\u\^ dx 

2 q Jrn 


< 


(/ _ 


1 


/ c(x)\u\‘^ dx = Iq. 

q Jrn 


□ 


Let us introduce 


ai = inf 


(/ - A)"/2^ 


> 0 . 


The main result of this section reads as follows. 
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Theorem 4.3. Under our assumptions on b and c, if 

(4.2) 

then equation (4-1) possesses a nontrivial solution. 

Proof. We follow the ideas developed in [9], By Ekeland’s variational principle, there exists a 
minimizing sequence {un}n C A4x such that lim„^+oo J\{un) = I\ and lim„_j._|_oo D,J\{un) = 0. It 
is readily seen that the sequence {un}n is bounded in so that we may assume without 

loss of generality that Un ^ u weakly in L“’^(]R'^). Let us set 

ttoo = lim limsup / \un{x)\^ dx 

R^+co n^+QO 7CB(0,R) 

Poo = lim limsup [ \{I — dx. 

R^+co n^+Qo JCR{0,R) 

We claim that 

(4.3) OliOolP < Poo 

(4.4) limsup / \un{x)f dx = / \u{x)f’dx + aoo 

n^+oo Jm.^ 

(4.5) limsup/ |(/— A)"/^ttnp dx > / |(/— dx +/3oo. 

n-j.+oo JrN 

Let us fix a smooth, positive cutoff function ip that equals one on a neighborhood of the origin. 
Define ipR = ip{-/R) and (pn = 1 — (pR. By definition of ai, we have 

Oll\\p>RUn\\l < ||(d - A){(pRUn)\\l 

By Lemmas 2.4 and 2.5 in the Appendix, 

(I - Ar/^i<pRUn) = (/ - Ar/\n -{I- Ar/^ipRUn) = 

(/ - A)“/\„ - ipRil - A^/^n + On{l) = (1 - ‘Pr){I - A^/^n + On(l) 

where 0^(1) 0 as n ^ +oo in L^(M^). Therefore 

aiWPRUnWl^ [ |(,5i?p|(/-A)"/^u„pdx + On(l) 

Letting n —>■ +oo and R —>■ +oo, we get (4.3). The relations (4.4) and (4.5) are easy and we omit 
their proofs. 

To complete the proof, we need to show that Uoo = 0. We argue by contradiction, and suppose 
that aoo > 0. Let us consider once more the function (pR. Since DJ{un){p>RUn) ^ 0 as n ^ +oo, 
we find that 

Poo P: ^botoQ. 


Combining with (4.3) yields 
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Moreover, 

J{Un) - ]rDJ{Un)Un = f ^ / h{x)\Unix)Y’ dx + f ^ A 

2 \2 pj Jrn \2 pj , 


> 


{P - 2)A 

2p 


c{x)\Un{x)\'^ dx 

N 

h{x)\un{x)Y‘(pR{x) dx, 


and letting first n —?■ +oo and then R —)■ +oo we have that 

p — 2 


h > 


2p 


-Xbar 


Together with (4.6) we get the contradiction 

p — 2 f ai 


h > 


2p \\h 




2p 


r' (a5) 


2 

2-p 


□ 


Remark 4.4. The lack of precise information about solutions to the limit equation 


{I — A)"n = b\u\^ ‘^u in 




and particularly the lack of a uniqueness result does not allow us to state the assumption (4.2) in 
terms of the ground state energy of the associated Euler functional. 

5. Concave-convex nonlinearities and sign-changing potentials 
In this section, following [5,20], we consider the equation 
(5.1) {I — A)°'u = \b{x)\u\^~'^u + c{x)\u\'^~‘^u in 

where N > 2a, A > 0 and and l<p<2<q< We assume that both |6| and jcj satisfy 

the compactness conditions (K). 

Remark 5.1. If —>■ 0 weakly in L“’^(R^), it follows from Proposition 3.5 that 


b{x)\un{x)\^ dx 


< 


/RJV 


\b{x)\\Un{x)\^ dx 0 


as n —>■ -|-oo along a subsequence. Therefore lim„_>,_|_oo /^jv h{x)\un{x)Y‘ dx = 0. By the same token, 
lim„^+oo /rjv c(x)|nn(x)|'' dx = 0. 


Definition 5.2. We say that a continuous function w changes sign if both sets {x S R" \ w{x) > 
0} and {x G R-^ | w{x) < 0} are nonempty. 


Let us recall from [15] that (5.1) is equivalent to the Neumann system 
— div (y^“^“Vu) -|- = 0 


m 


oN 


I —y^ = Ka {Xb{x)\u{x,0)\^ ‘^u{x,0) + c{x)\u{x,0)\'^ ^u(x,0)} on R^ x {0} 

Weak solutions to (5.2) correspond to critical points of the Euler functional Ex ■ 


dehned by 
Ex{v) = I 




(^\Vv{x,y)\^+ v{x,yf^ dxdy 


A 

p ArJv 


6(x)l trv{x)\^ dx — 


1 

q Jrn 


c(x)j tr v{x)\'^ dx. 
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where is the completion of with respect to the weighted norm 

(|Vi;(x,y)p+ i;(x,y)^) dxdy 

and tr: ^ is the continuous trace operator defined in [15, Proposition 

6.2], In particular 

(5.3) Kail tri;||^„,2(KAf) < II^II^1 (r^+1) ^ ^ 

Remark 5.3. Up to a rescaling, we can (and will) assume that Ha = 1. 

The Nehari manifold associated to (5.2) is 

A6. = {u E H\R^+\y^-^^)\{0} I DEx{v)v = o} 

or the set of those v E , y^~‘^°‘) \ {0} satisfying 

(\Vv{x,y)\'^ + v{x,y)‘^") dx dy = X 6(x)| tr i;(x)|^ dx + / c(x)| tr i;(x)|'^ dx. 

^ ^ jrk 7rk 

To each v E id^(R+^^, i/^“^“) we attach its fiber map (fy : [0, +oo) ^ M defined by (pv{t) = E\(tv). 
By a direct calculation, 



\h^ 


\yl-2o.) 


, 1 - 


,Af+l 


Dipy{t)=t f|Vx(x, y)p + x(x, y)^^ dxdy — XE ^ f 6(x)| tr x(x)|^ dx 

— [ c(x)| tr x(x)|'^ dx 

Jr^ 

and 

D‘^(py{t)= f^^y^~‘^‘^(^\Vv{x,y)f + v{x,yf) dx dy - {p - l)XtP~‘^ f b{x)\tTv{x)\P dx 

— (r — l)t'^“^ / c(x)| tr x(x)|'^ dx 

JRK 

Remark 5.4. Clearly, tv E A/a if and only if Dipy{t) = 0. In particular, v E A/a if and only if 

D(py{l) = 0 . 

We decompose 

WA=Ar,^uA/^uW+, 


= {u E Wa I T>V^(1) < o} 
Wa° = {ueWa I IlV.(l) = o} 

= {u E Wa I T>V^(1) > o} . 


where 
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Let US define 


G ^") | J c(x)| tr i;(x)|'^ dx > o| 

C- = L e \ [ c(x)|tri;(x)|''(ix < ol 

I Jrn J 

= IV G I [ c(x)| tr i;(x)|'^ dx = ol 

I JRJV J 

= |x G I J 6(x)| trx(x)|^(ix > o| 

B~ = Iv G \ [ 6(x)| tri;(x)|^(ix < ol 

I Jrn j 

= Iv G H\R^+\y^-^‘^) I [ 6(x)|tri;(x)|P(ix = ol . 

I JRff j 

Lemma 5.5. {Ij If v G B~ n C~, then no multiple of v belongs to M\. 

(2) If either v G B'^ n C~ or v G B~ n C~^, then there exists one and only one t{v) > 0 sueh 
that t{v)v G J\f\. 

Proof. Both statements can be proved by an elementary inspection of the fibering map ip^- □ 

The case v G B'^ n C~^ requires some additional care. 

Lemma 5.6. There exists Aq > 0 such that A < Aq implies (/?„ > 0 for all v G , y^~‘^°‘). 

//A < Ao and v G B^ n C~^, then possesses exaetly two eritieal points. 

Proof. Let v G H^{R^~^^,y^~‘^^) be such that /j^jv c(x)| tr x(x)|^ dx > 0. We introduce 

F'v{t) = pr [ (\Vv{x,y)f + v{x,yf) dxdy-— [ c(x)| tr x(x)|‘? dx. 

2 JR^+i ^ ' Q Jr^ 


C- = G 


= \vGH 


lon).^+l 2 q\ 


B+ = G 


B- = <!x G 


B^ = Iv gH 




Since 


DFJt) = 


(|Vx(x,y)p+ x(x,y)^) dxdy-t*? ^ [ c(x)| tr i;(x)|''dx, 

JR^+l ^ 2 J^N 


the function K, attains its maximum at 


J^N+iy^ {\V{x,y)\^ + v{x,yf) dxdy^- 
/jjAT c(x)|trx(x)|9dx 


Moreover 


^\ (■fR^+^y (|V(x,y)|2+ i;(x,y)2) dxdy , 

F^{t*)=(---] ^5 - 1 

V2 qj y c(x)|trx(x)|‘?dx)^ 

D‘^Fy{t*) = {l-q) / y^"^" (|Vx(x,y)p + x(x,y)2) dxdy < 0. 

JM. I 


,2 q 


Urn c{x) \ tr V{x)\<i dxY 


Let Sg be the best constant for the inequality 


|trx(x)|‘?dx) " < yi (|Vx(x,y)|2 + x(x,y)^) dxdy 
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that follows from (5.3) and Theorem 2.2. Then 


with 6 independent of v. Moreover, 
{t*)P 


1 1 


1 

q-2 


2 qj \\\b+\\o.S^^ 


= 5>0 


p JM.N 


< 


b{x) \ tr v{x)\P dx 


-SP 


'J^N+iy^ (|Vi;(x,?/)p + i;(x,y)^) dxdy\<t 
/gjv c(x)| tr v{x)\^ dx 


P 

-2 




P 


> qp 
Dp 


(|Vx(x,y)p+ x(x,y) 2 ) dxdy) ^ 


‘l\ 2 (,- 2 ) 


Urn c(x)|trx(x)|'?dx)" 


/ 


i QP 

Jp 


p 

2q \ 2 


q-2 


F,{t *)-2 = cF,{t*) 2 . 


□ 


Therefore 

Pv{t*) > Fyit*) - XcFy{t*)2 = Fy{t*)2 (F„(r)^ - Ac) > 62 - Ac) . 

1 2—p 

We complete the proof by choosing A < c “^5 2 = Aq. 

Corollary 5.7. Let A < Aq. There exists 5i > 0 such that E\{v) > 61 for every v S -Xf\ ■ 

Proof. Indeed, if x G ^ then cpy has a positive global maximum at t = 1, and in addition 

/ 6 (x)| trx(x)|^dx > 0 . 

Jrn 

Hence 

Exiv) = Ml) = Mt*) > EM)"^ {Ey{t*)"-^ - Ac) > (( 5 ^ - Ac) > 0 

provided that A < Aq. □ 

Corollary 5.8. If X & (0, Aq), then = 0. 

Lemma 5.9. Let v be a local minimum point of Ex on either or Aff) such that v ^ Af^- Then 
DEx{v) = 0. 

Proof. By the Lagrange Multiplier Rule, DEx{v) = pDIx{v) for some // € M, where we have set 
Ix{v) = DEx{v)v. Hence DEx{v)v = pDIx{v)v = pD'^ipy{l) = 0. Since v ^ Af^ we must have 

p = 0. □ 


Lemma 5.10. The restriction of Ex to Af is bounded from below and coercive. 

Proof. Indeed, for any v € Afx, we have 

E,(V) = (I - I) = 

and we conclude since p <2. □ 
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Lemma 5.11. If X < Aq, then Ex attains its minimum onM^ ■ 

Proof. We k n ow from Lemma 5.10 that Ex is bounded from below on Therefore there exists 
a sequence {vk}k C with 

lim Ex{vk) = inf Ex > —oo. 

k^+QO _xf+ 


Again from Lemma 5.10 the sequence {vk}k is bounded. We are then allowed to suppose, without 
loss of generality, that it converges weakly to some vx in Let us pick v S 

such that fj^N b{x)\ ticv{x)\P dx > 0; then by Lemma 5.5 there exists ti = ti{v) > 
0 such that tiv E and Ex{tiv) < 0. This implies that inf_^+ Ex < 0. On A/a we have 


and thus 




By Proposition 3.5, Remark 5.1 and a semicontinuity argument, we let k +oo and get 


We finally claim that Vk vx strongly. If not, then 

J^N+i (|Vu(x,y)|^ + u(x,y)2) dxdy < limh£ (|Vufc(x, y)p + Ufc(x, y)^) dxdy. 


Let us recall that 


=t (|Vufc(x,y)p+ Ufe(x,y)2) dxdy 

— Xt^~^ [ b{x)\tTVk{x)\^ dx — [ c(x)| tr Ufc(x)|'^ dx 

JK-v Jrjv 

and 

DipvA'^)=t 2 /^“^" (|V'f^A(x,y)|^+ UA(x,y)^) dxdy 

— [ b{x)\tii:vx{x)\^ dx — [ c(x)| tr ua(x)|'^ dx. 

Jrn Jrjv 

Let tx = tx{vx) > 0 be chosen in order that txvx E this is possible by Lemma 5.5. For k ^ 1, 
we have Dipy^{tx) > 0, which yields tx > 1. But then, by definition of tx., 

Ex{t\vx) < Exivx) < , lim Exivk) = infEx- 

k^+oc xf+ 

This is impossible, and we must have Vk vx strongly in as k ^ + 00 . Since 

Mx = 0 , there results vx E and the proof is complete. □ 

Lemma 5.12. If X < Aq, then Ex attains its minimum on Nff. 
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Proof. We know from Corollary 5.7 that Ex > 0. Let {vk}k C Mff be a minimizin sequence, 

i.e. 

lim Ex{vk) = inf E\ > 0. 

fc—>-+CXD 

By coercivity on Mff we may assume without loss of generality that {vkfk converges weakly to 
some vx in . Now, 

Exivk) = (^ - ^) + (^ - J) J^^c{x)\tTVk{x)\Ux. 

By the assumptions on c. Proposition 3.5 and Remark 5.1 

lim / c(x)| tr Ufc(x)|'^ dx = / c(x)| tr UA(a;)|'^ dx. 

k^+ooJ^N JrJV 

Hence /j^Af c(x)| tr ua( 3:)|'^ dx >0. As a consequence by Lemma 5.5, the fiber function possesses 
a global maximum at some t = t{vi) such that tvx S -^x ■ other hand, Vk € Mx already 

implies that 1 is a global maximum point for namely p>Vf,{t) < for every t > 0. 

If {vk}k does not converge strongly to vx in then 


E- 


■ itv\) = ^ [ ^"(|VuAp + ui) dxdy - — j 6(x)| tr i;a(x)|p dx 

c(x)| tr vx{x)\‘^ dx 


E 
q Jr^ 


P f 

< lim inf — / 

fc^+oo 2 jRf+^ 


(^\S7vkf + vl'^ dxdy - 6(x)|trufc(x)|Pdx 

E f 

-/ c(x)| tr Ufc(x)|'^ dx 

q JR^ 


< lim Ex{ivk) < , lim Ex{vk) = inf Ex, 

K —^“|~00 K —^“|~00 yy /” 


a contradiction. We have proved that Vk vx strongly in as k ^ +oo, and 

Afx = 0 hnally implies that vx G ■ □ 


We are ready to prove our main result. 


Theorem 5.13. Let N > 2a and l<p<2<q< 2N/{N — 2a). Suppose that both |6| and 
|c| satisfy the compactness conditions (K) and that b and c change sign according to Definition 
5. 2. Then there exists Aq > 0 such that for every 0 < A < Aq equation (5.1) has at least two 
non-negative solutions. 


Proof. We only have to prove that the two minima of Ex on J\f)( and on respectively can be 
chosen to be non-negative. To this aim, we let 


f{x,t) 

f+{x,t) 

E+{x,t) 


6(x)|t|^ -|-c(x)|t|'^ ‘^t 

f/(x,t) ift>0 
1 0 if t < 0 



/+(x,s)ds. 
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To the functional 

I^^F+ix,vix))dx 

we can apply all the previous lemma, and thus for every 0 < A < Aq the functional possesses 
two critical points ua,i G and v \^2 G -^x ■ Now, denoting by the negative part of ua,i, we 
have 

0 = (V0A,1 ■ dxdy 

Therefore ua,i = > 0. Since the same computation holds with vx ,2 in place of vx^i, the proof 

is complete. □ 


6. Perspectives and open problems 


As we have seen, the Bessel operator (/ — A)" shares many features with the standard fractional 
Laplacian (—A)". However, the most striking difference between these two operators is that the 
latter has some scaling properties that the former does not have. As should be clear from (2.1), 
(2.2) and (2.3), the Bessel operator is not compatible with the semigroup M+ acting on functions 
as X i-A u(s“^x) for s > 0. In simpler words, the Bessel operator does not scale. 

In our opinion, a very challenging problem is that of finding solutions to the fractional scalar 
field equation 

(6.1) {I — A)°‘u = g{u) in 

under the so-called Berestycki-Lions assumptions on g. In [3, 4] the case a = 1 was studied 
under very mild assumptions on the nonlinear term g. In particular, no Ambrosetti-Rabinowitz 
assumption must be imposed, and no monotonicity assumption like in (f3). 

However, the main tool used in [3, 4] consists in a clever exploitation of the semigroup action 
s-ku: X i-A u{s~^x) for s > 0. Indeed, solutions are constructed (roughly speaking) by solving the 
minimization problem 

inf I [ |Vup + u‘^ \ [ G{u) = 1 

[Jm.n Jrjv 

where G is the antiderivative of g. Then a rescaling from u to a suitable s-ku produces a solution of 

(6.1) by absorbing the Lagrange multiplier. Let us propose a closely related question by considering 
the special case of equation (1.1) 

(6.2) {I - A)^u = \u\y-‘^u inR^. 


We continue to assume that 1 < p < 2* . In [17] it is shown that (6.2) has infinitely many solutions. 
We claim that it possesses a radially symmetric ground state. Indeed, we consider the quotient 


5 = inf 


(l - A)“/\) 


dx I u G L 


0 , 2 / 




|u(x)|^ dx = 1 


To overcome the lack of compactness in R'^ we can work in the subspace L[[![^(R'^) consisting 
of radially symmetric elements of L"’^(R'^). The fact that S is attained is now an immediate 
consequence of the compact embedding of L]]![j(R^) into LP(R^), see [32]. To prove that the 
minimizer of S can be chosen to be radially symmetric, we use the following inequality: it is 
probably known, but we could not find a precise recerence in the literature. 
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Proposition 6.1. Assume that u S L“’^(]R^), and let u* the usual symmetric-decreasing re¬ 
arrangement of u. Then 


(6.3) 



(I - u* 


2 

dx < 


{I 


dx 


Proof. We adapt the proof given by [25, Appendix A], see also [24]. By [15, Eq. (7.5)], 
Ka [ (-A + l)"/^rt(x) 

JRN 


dx = —2a 


i?(vWTTt) -1 


t 


\Fuicrd(, 


where Ka is an explicit positive constant and i? solves the ordinary differential equation 


1 — 20 

r +- d' - 1 ? = 0, i?(o) = 1. 


Moreover, 


(6.4) f dC = t f (u{x)P{t,-)*u{x)-u{xf] dx, 

JrN t^» J^N \ / 


where 

_ Ai^ + 2Q: / _\ 

P{t,x) = C'f^at‘^°‘ + K n+^2c K/|xP + , 

P{t, •) * u{x) = J^N P(t, X — y)u{y) dy and (7)^- ^ is another explicit positive constant. Since KN+ 2 a 

is positive and non-increasing (because Ku{s) > 0 and Kl{s) = —^Ku{s) — Ki,-\{s) < 0 for all 
s > 0), we can apply an inequality by Riesz [24, Eq. (3.9)] and conclude that the right-hand side 
of (6.4) decreases if we replace u by u*. We achieve (6.3) by letting t —>■ 0. □ 


The proof of the following result is now straightforward. 


Theorem 6.2. Let 1 < p < 2* . Then eguation (6.2) possesses a radially symmetrie ground state. 

Remark 6.3. By the results of [27], if u = G 2 a*u^ with u E and r > max {p, N{p — 1)/(2a)} 

then u is radially symmetric and decreasing about some point. The previous result is weaker but 
much easier to prove. 


What happens if we replace the power in (6.2) with a more general nonlinear term 

f{u) with subcritical growth at infinity? It is not hard to check that we can solve the minimum 
problem 


5 = inf 


/rJv 


(/ - A^/'^u) 


dx 


E L 


a,2tJa,N\ 


/RJV 


F{u{x)) dx = 1 


in the framework of radially symmetric functions, but getting rid of the associated Lagrange 
multiplier is a complicated task. This difficulty is again caused by the lack of scaling invariance 
for the Bessel fractional operator. 

Eor very similar reasons, it seems that the Bessel operator (/ — A)“ does not produce strong 
variational identities like the local Laplacian or the fractional Laplacian. Since identities like 
Pohozaev’s are a consequence of the action of the very same semigroup s -k u, we can imagine 
that some troubles arise. Indeed, Felmer et al. proved in [17, Proposition 5.1] that the following 
pointwise identity holds for every ip E 5(M'^): 


(6.5) (/ - A)“ (x • V<p) = X • V [(/ - A)“<p] + 2a{I - A)“(p - 2a{I - A)“-V- 


Here the crucial remark is that such a formula for an “integration by parts” involves the inverse 
Bessel operator (/ — A)““^, since a — 1 < 0. In other words, a direct attempt to extend the 
Pohozaev identity leads to terms that are different in nature from those appearing in the problem. 
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Indeed, if (I — A)“ii = f{u) and u decays sufficiently fast at infinity, we can justify the following 
formal computation, see [17] for a similar result: 

[ {x \ Vu){I — A)^udx = [ {x \ Vu) f {u{x)) dx. 

Setting 

(/) = [ (x I Vu)(/— A)"u(ix, {II) = f {x \ Vu) f {u{x)) dx, 

we deduce that 

{!) = {2a — N) [ u{x) f {u{x)) dx + N [ F{u{x))dx — 2a [ u{x){I — /S.)°'~^u{x) dx 

JRJV jRAf 

and 

{II) = -N f F{u{x)) dx. 

Jrn 

Therefore the solution u satisfies the variational identity 

(6.6) 2a [ u{x){I —/S.)°'~^u{x) dx = 2N [ F{u{x)) dx — {N — 2a) [ u{x)f{u{x))dx. 

The left-hand side of (6.6) is finite if and only if u E L““^’^(M^), and this is certainly true since 
0 < a < 1 and therefore L“’^(M'^) C L"“^’^(M^), see [26, Equation (7.3)]. Anyway, the identity 
(6.6) is structurally different than the usual Pohozaev identity in the local case a = 1. 

We observe that for the standard fractional Laplacian (—A)“ the analogous of (6.5) reads 

(—A)“(x 1 Vu) = 2a{—A)^u -|- {x \ V(—A)“u), 

see [31, page 602], so that we get a more useful variational identity. 

It would be interesting to understand if classical results like those proved in [3,4] using varia¬ 
tional identities can be somehow extended to Bessel fractional operators. 
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